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ABSTRACT 
Considering the absorption of terrestrial radiation by carbon dioxide 
and water vapor, coefficients of eddy conductivity are computed Гог heights 
in the range 4 to 320 meters from data collected during the Great Plains 
Field Program (August and September 1953.) Power-law profiles with height 
are fitted ۳ the computed coefficients and values obtained from the fitted 


profiles correlated with a stability parameter. 
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1. Introduction. 

An important step forward in the theory of the turbulent transfer of 
properties in the lower layers of the atmosphere was taken when the con- 
cept of Austauschkoeffizienten or "exchange coefficients" was proposed 
by Wilhelm Schmidt in 1925. Since that time several attempts have been 
made by бейе workers to measure these coefficients through experi- 
mental observations in order to test the theory of turbulent exchange. 

In 1950 Cooley yes. / made determinations of eddy conductivity near the 
earth's surface from data gathered at Manor, Texas. About the same 

time Jehn and Gerhardt /9 / made measurements of eddy heat flux up to 
about 120 meters at Manor and Gila Bend, Arizona. In addition to these 
were the observations of Cramer and Record /6 / near the earth's surface 
in 1553, and a small number of observations due to other workers in the 
period 1955-1957. These observations have been summarized by Priestly 
/iV, end no attempt will be made here to review all of the results. 

E 

The primary reason that the number of such observations has been 
small is that the accurate determination of micrometeorological parameters 
is exceedingly difficult and very expensive. For this reason the quality 
and volume of such data obtained in the past have been far below the 
level of desirability. 

The publication of the "Proceedings of the Great Plains Turbulence 
Fieid Program, 1 August to 8 September 1953, O'Neill, Nebraska" made 


new and better data available for the study of the turbulent transport of 


properties in the atmosphere's lowest layers. In order to obtain new 





evidence on the nature of the turbulent heat and momentum transports 
and their relationships to one another, it was originally decided to 
initiate thesis projects at the United States Naval Postgraduate School. 
These projects to be undertaken were to determine the coefficients of 
(1) eddy conductivity, Кн, (2) eddy viscosity, KM , and (3) the ratio 

oi the two alte ; K Ku from this new micrometeorological data. 
In particular this present study involves the determination of the coeffi- 


cient of eddy conductivity and the examination of its variation with 


height and stability from the values determined. 





дн The Great Plains data. 
The data collection program for the Great Plains project was broren 
iown into seven time-observation periods with a multitude of parameters 


being observed at nearly equally spaced observation times within each 


For this study the seventh observation period (7 and 8 September) 
w35 used end was considered to be the most suitable of the periods on 
the basis of a combination of two factors; (1) completeness of data uti- 


lized in the study and (2) a minimum of advection during the period. As 


c 


an aid in the interpretation of the results the following quotation / 10/ 
is provided concerning the synoptic conditions in the area during the 
seventh period: 


On 7 September , the high-pressure area east of Nebraska moved 
eastward while a trough of low pressure approached from the west. The 
winds were southerly and initially light, becoming 15 to 20 mph. The 
Skies were generally clear over Nebraska with some scattered cirrus 
during the morning of the 8th. Cumulonimbus clouds and lightning were 
visible to the north of the test site from the afternoon of the 7th, until 
the early morning of the 8th. The seventh and last general observation 
period began at 0405 CST on 7 September and continued until 1305 CST 
on ë September under relatively favorable conditions. 


Of the many tables of data included in V, the following summary 


snows those which were used in this study: 


Га) 


Table number 1 | Data used 
1.1а 7 September, 1900, CST, West and 


South components of the geostrophic 
wind at indicated pressure levels. 
Values estimated from local synoptic 
charts. 


2 


Table where numbers of the form 1.1 etc., refer to tables in /10/ 
^W ere times are omitted, data were available for each soundinc d" 
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“% 
Table Number! Data Used “ 


172 Station pressure and reletuee Nu halty 
o.là Soil temperature at 0.5 em below the 


surface of the ground (taken as being 
representative of the ground tempera - 
ture .) 


4.2 * Mean air temperature from aspirated 
thermocouples (MIT, 15 мк” mesan 
at 2.0, 4.0, 6.6 and 15.0 a above 
the ground, UCLA, average Six 
readings taken every 3 minutes at 
various levels up to @.0 meter 
shielded thermocouples СН, cs 
of 20 readings at each levei during a 
5-minute interval et variouós levels up 
to 6.4 meters) and shielded thermisters 
(UT, 10-second mean: at variou. levels 
up to .8 meters.) 


5.28 Mean cross product ot vertical кы 
velocity component and gir tempera- 
ture multiplied by density and cias 
heat of air for 3.0, 6.0 and 12.0 meters 
above the ground, value: missing for 
1235 and 1435, 


6.4 Radiosonde pressure, temperature and 
relative humidity at signiiicant levels. 

G.S Air temperature and mixing ratio ae alti- 
tudes from aerograph (Len) data. ше 
at start, top level and end of Might is 
indicated, 

7.48 Sensible heat transfer to tle air from 


the earth-air interface sccording to 
theoretical models o: Suomi and Lettau. 


During the observation period, data were recorded bi-nourly tor even 
hours of average local time (CST plus 35 minutes). Sufficient time and 
Jata were available to make Ku computations for six soundings, 0535, 0855, 


1232, 1235, 1435 and 1635 on 7 September. Since no night-time зегостари 


observations were made, only daytime determinations were possible. 
4 





2.1. Instrumental errors. 
Contributors to the Great Plains project have given the following 
values for instrumental errors inherent in the data: 


т о 0 3 г 4. 
l. Micrometeorological absolute temperature errors (Table 4.2) 


«> © 


a) 0.2°C (UT) 
sb) 0.1°C (MIT and UCLA) 
c) 0.05°C (JH) 


2. The L-20 aerograph configuration? temperature and relative 
humidity errors. 


a) Temperature error: 
(1) 1°C absolute, for any level. 
(2) Relative error between levels, 0.19c. 


b) Relative humidity error: 
(1) to 10% absolute. 
(2) Relative error between levels, 3%. 


In 29 cases a comparison of MIT temperature values at 16.0 meters 
to the L-20 data at 17 meters gave a mean difference, MIT minus L-20 of 
- 0.10 + 0.119С. 

For a complete picture of the instrumentation, analysis of instru- 


mental errors and a description of the test-site the reader is referred to 


lapor an interesting comparison of observed micrometeorological tempera- 
к see fig. 1. 
The L-20 configuration utilizes the KS-4 Kollsman aerograph, which is 
an automatic device for the automatic recording of temperature and rela- 
tive humidity. 
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За Background, methods and equations. 


The coefficient of eddy conductivity, Кн, is defined by: 


_ EN (1) 


Cs c 
CGE 


aye 


where © is the turbulent heat transfer in the vertical. It may be observed 


from this definition that K., is always a positive quantity. When 98 is 


H 
negative (unstable conditions), Q will be positive and directed upward. 
Conversely for stable conditions (92 е о). Q will be negative and directed 
downward. 
In this study starting values of Ка at 4 and 8 meters ere computed 
bv using eq. (1) with Q given by the average value in the laver 2-12 
meters. For values above the starting level the following method is utilized. 
The horizontal turbulent exchange of heat is usually very small 
compared to that in the vertical, and vertical advection is usually con- 
sidered to be negligible in comparison to horizontal advection. Tne fact 
that clouds did not occur in the layer under study rules out a latent neat 
contribution, and the small height-variation of density can be overlooked 


for the layers to be considered. 


Then the equation of turbulent diffusion for heat becomes: 


(98) .- v8 + 2. («9€ (2) 
01-7 8 
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Now returning to eq. (2) and substituting eq. (3) we get: 
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3.1. The water vapor problem. 

The water vapor radiation problem was approached using a method 
presented by D. L. Brooks ya / and based upon an adaptation of the 
flux divergence equation of Bruinenberg: 


۱ i Ша | vo gu 
T | ki Ki p ЈИ | > 人 вт“ Es; 
0 h m n ۸ 


\ot/wi Ср \ ۵ 


| а Y 
where qis the specific humidity at the reference level at which e 


Jl 
۱ 2, | 
is to be computed; (Es) 1s the Elsasser pressure correction factor; 

A T, | ۱ 
and ЛО ' jig the change in the Stefan-Boltzmann black body flux between 
successive levels of temperature | in the sounding. The summation 
index n denotes successively higher levels in the atmosphere above the 


reference level and the index m the layers below the reference level. 


mS 


^e partial derivative dis the slope of the curve of emissivity of water 
Qi | 


= 
Jt. 


CES 4 2 Є. 
vapor, plotted against corrected optical depth and the term ۱ У is 


evaluated at the top of the sounding. 

The sctors 96 are evaluated as finite difference slopes from 
Brooks' curve! . In addition to the curve Brooks / ШЕ also presents 
tabular values of the mean slope of the emissivity curve for small incre- 
ments of optical path-length. 

If the convention in using eq. (7) is that des is always positive 


aw 


when temperature decreases with height above the reference level, cooling 


^*Not shown here but presented in / 1/. 
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at the reference level results from the layer contributions above. 
Similarly when the temperature from the ground decreases up to the ref- 
— level, the contributions from below yield warming at the reference 
level. The second term on the right, representing radiation divergence 
from the top oi the sounding, gives a cooling contribution at the reference 
level. The actual conventions used in working with the tables are some- 
what different, however, and were designed for maximum efficiency in 
machine computation. 

In the application of the method, the initial step is the computation 


of the corrected optical depth, w, given by” 
| p. las Ө 
af o Г. 


were p. and p are the layer boundary pressures and p, is a standard 
O 2 5 


pressure, 1000 mb. 


ts 


‘g, (8) is approximated for a given layer by using mean values of 


qg and p, and a finite pressure increment Др: 
əй (p N (9) 
la wW = 24 АР „Е 


When increments of optical path length, Aw, are calculated tor 
each layer, values of 95; may be obtained by interpolation in the table, 
using values w and w + Àw in entering the table. 

When this work was begun the decision was made to perform the 
computations by hand methods rather than by electronic digital computer. 
This decision was in part based on the fact that use of the then available 
computer, the National Cash Register Company Model 102A, would have 


10 





necessitated lengthy computation times exclusive of the time required 

for programming. This observation is based on the fact that there are 
over 1000 entries in Brooks' table and only about as many storage cells 
in tre computer. If the 102A had been used it would have been necessary 
to store the entire table on magnetic tape, which would virtually have 
eliminated the computation time-advantage of the computer, This diffi- 
culty may now have been overcome with the advent of the Control Data 
Corporation model 1604 computer. The larger storage capacity (over 
32,099 cells) and extremely high speed of this new computer would lend 


itself well to Brooks' tabular method. 
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3.2. The “arbon “ioxide y adiation problem. 

For the carbon dioxide part of the radiation problem it was decided 
to employ a simple set cf working equations previously used by Brooks 
ГУ; in stratospheric investigations of carbon dioxide radiative transfer. 

The downward flux of infrared radiation at a level in the atmosphere 
and for a given frequency band (12.2 - 17.5,4 for CO )) is given by the 


Schwarzschild equation integrated over all angles: 


И Bo: 
FY) = FLO) (om) +) Bb AC. че 


956۲۵ ۲۱۸1۵ the downward flux of diffuse radiation at the reference 
level, where w=w, ; МА is the optical path length of CO, above the level; 
F ¿ (0) is the flux down from the top of the atmosphere; E. (0,1) is the 


transmissivity of diffuse radiation by a slab of optical path w,; and 


1 
B, = „185 < т4 is the total black body Нива radiated by CO» at the approp- 
riate mean temperature of the үш layer. 
Ihe summation index i runs through n layers from the top of the atmosphere 
to the level in question. AT; is the increment to the flux transmissivity 
at the level Wı contributed by each layer. 
Similarly the upward flux F ^ (w) is: 
FIM) = Fale) (mA, A) + а B. A Се; (10) 


where F Á (м) is the upward flux at the ground; and Da (w, Wo) is the 


transmissivity (assumed independent of temperature) of the slab of 


I-A carbon dioxide contribution equal to 18.596 of the black body rad1ation 


is quoted by various authors. 
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optical thickness wg-w,, or from the ground to the reference level. The 
summation j-1 through m runs from the ground to the level W]. 

In using these equations it is assumed that each layer can be app- 
roximated by an isothermal layer whose absolute temperature is the actual 
mean fourth rọot of the temperature of the layer raised to the 4th power. 
Since the soundings do not run to the top of the atmosphere, F í (0) was 
computed at the top of each sounding. 


The net flux at the reference level (м1) is then simply: 


By ie) = Grp) ea) aa 


and the flux divergence, AF, , is taken as being the difference of net 
flux at two adjacent levels, the assumption being that the net flux 
divergence is uniform between the two layers. The hourly-temperature 


change in terms of the flux divergence in the layer is then: 


ОТ) талады (13) 
(91 CO2 l ^p 





where flux is in EEN min and pressure is in millibars. 

For the determination of flux-transmissivity 7;., several choices 
were available as to method. СаПепдет /4 /Gives formulas for T based on 
a breakdown of the 12.5 - 17.5 band into five smaller wave length 
intervals. Burch and Howard PER and Howard, Burch and Williams 
/8 / give two empirical formulas based on laboratory data for the total 
absorption of the 12.2 - 17.5 44 band: 


"Strong" fit; "m dy: 365+ 55 Хед. № + 4T Leg (Io if) (14) 
"Weak" fit: SA, dp= 3-16 w 2 (Рев) 4 (15) 
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Б Ado ГАРА 
It follows then that A = 777; “Rese am! 


and Z = 7 一 4 

In (14), (15) w is absorber concentration in atmosphere-centimeters; p 
is atmospheric pressure in mm Hg; and ó is the partial pressure of CO, 
(assuming a constant ratio by weight, .00046, of the atmosphere.) 

The Howard et al expressions (14) are for the "strong" fit or rela- 
tively large values of w and (15) for "weak" fit or relatively small values 
or w. The cut-off value given is 20% absorption, the log expression 
applying for A> 20% and the exponential expression applying for < 20%. 

Both of the methods discussed above have the disadvantage of 
extremely lengthy computation-time for the data under consideration unless 
an electronic digital computer is used. Furthermore, the empirical equa- 
tion of Howard et al has the disadvantage of a discontinuity between the 
two equations near absorptivities of 20%. A third and faster choice pre- 
sented itself in the form of a curve giving fractional absorption A vs pres- 
sure x path length (mb-cm-atmos) given by Brooks ОБ, апа credited by 
Brooks to Burch and Howard /3/ (see fig. 2.) This curve had the advan- 
tage of affording greater facility in conjunction with hand computations , 
and it had no obvious discontinuities. A difficulty presented itself, 
however, after the CO, computations were well underway when it was 
discovered that use of the Brooks curve yielded positive CO, radiative 
flux divergences near the ground where large negative water-vapor flux 
divergences were obtained. Since this occurred near the level where 


absorptivities of 20% might be expected, the difficulty was overcome by 
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changi ng from the curve to the Howard et al "weak" fit expression for 
curve absorptivities less than 20%. 

Although Brooks credits the curve to Burch and Howard Гай it is 
not entirely clear how he obtained his curve. The Burch and Howard curve 
is similar to that of Brooks except that their ordinate is percent trans- 
mission and the abscissa parameter is absorber concentration (cm-atmos). 
The Burch and Howard curve was obtained from data based on the use of 
the two empirical formulas, the discontinuity being smoothed out in the 
construction of the curve. Thus, regardless of how the smoothing was 
done in obtaining Brooks' curve it is clear that it is based on the empirical 
expressions of Burch and Howard. On this basis, changing from the curve, 
fig. 2, to the "weak" fit formula where necessary is justified, the discon- 
tinuity being accepted. 

For the sake of comparison, values of percent absorption due to the 
"weak" fit formula are plotted in fig. 2 along with Brooks' curve. In the 
latter, a value of the "strong" fit formula corresponding to an curve absorp- 
tion of 2096 is shown. It appears that the smoothing in the region of 
10-20% absorption favored the "strong" fit curve. 

Since the fractional concentration of CO, by weight in the atmos- 
phere is nearly constant with height, the amount near any pressure level 
over a given vertical distance may be expressed in terms of a pressure 
difference. Thus every 4.2 mb of atmospheric air is equivalent to 1 
cm-atmos of carbon dioxide (ie, a layer of CO» of thickness 1 cm if 


reduced to S. T.P.). Taking the diffuse radiation through a slab to be 
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1.67 times that through a parallel radiating column, the equivalent 


Optical path is: 


E ied 3 о (сло ATMO 
= 247 fe "dp = 394 Др Em- arnes) (16) 


This is the absorber-concentration in the Burch and Howard expression. 


For use with Brooks’ curve, w is linearly pressure-weighted giving w as 


follows: 
| P» 
/ 0 (067 | = 
W = 244 ы pap = . 394 P др (17) 


with w' now in mb-cm-atmos and p the mean pressure of the layer. 
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4. Computations. 
4.1. The water vapor radiative temperature change. 

The radiative temperature changes due to water vapor were com- 
puted by use of the method discussed in section 3.1. Required are the 
knowledge of temperature, pressure, pressure differences, and relative 
humidity at the given reference levels, the reference levels being taken 
as 8, 17, 35, 100, 165, and 240 meters. Here, conditions at 16.0 meters 
are assumed to apply at 17 meters, so that in some cases, in order to 
smooth out artificial inversions in the temperature data between 16.0 
and 17 meters, it was necessary to shift the MIT temperatures so that 
the 16.0 meter MIT temperature would correspond to the 17 meter L-20 
temperature. 

The pressure, and layer pressure differences were computed by 
use of the hypsometric equation together with the known temperature 
profile. Values of specific humidity were obtained from the aerograph 
humidity data, radiosonde data, and station humidity. Unfortunately 
instrumental failure left a hole in the micrometeorological humidity data 
between the station values and the aerograph soundings. For levels 
between station values and aerograph soundings, linear interpolation 
with height was assumed valid in arriving at values of д for the reference 
levels. As noted, radiosonde data was used to extend the data above 
the maximum level of L-20 flight. This was done by "hooking" on to the 
nearest radiosonde reporting level after the highest reported L-20 data 


were exhausted. 
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The computed water vapor radiative temperature changes are 
given in table 1. Note the strong radiative warming tendency in the 
lower layers at midday , tending to change to cooling in the neighborhood 
of 100 m. This effect has already been noted by Cooley /5 7. 
Table 1. Infrared radiative temperature change (VC /Hr.) 
due to water vapor (times in CST). 


Reference Times 
level (meters) 0635 0835 1035 1235 1435 1635 


8 -0 13 1.25 2.00 3.05 223 1۳ ۱ 
Pp 20:07 Шы ay [57 0 92 0.10 
35 00 0.14 ШІ 0.13 ШЕ 510 0207 
oul 01409 0.09 0.07 =0 0.10 -0.14 
100 -0.14 0.03 0.02 022 оо -0 Ша 
165 =) al 0.09 x0 93 0 uU ie 
240 ЕО "TT -0.05 -0 05 -0.06 -0.05 
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4.2. The radiative temperature change due to carbon dioxide. 
ine са radiative temperature changes were computed by the 


method of section 3.2 and are shown in table 2. Since the mean black 


adm m 


4 4 
body flux oT for layers was needed in the CO, method, plots of of 


2 
versus height were made for the layers from 0 to 16 meters, and 0.185x 


ST ‘approximated by equal-area averaging. Above 16 meters, 0.1856 т“ 
was simply averaged between levels. 


Since the method yields Q9! 


А іп layers between reference levels, 


OF 
reference-level values were interpolated. On the average, co, radiative 
temperature changes were found to be approximately 23% of the water 

vapor radiative changes, a result which also confirms that of Cooley VAs, / 


in these layers. Note also the tendency for the radiative temperature change 


to change sign near 240 meters in a manner similar to that of water vapor. 


Table 2. Infrared radiative temperature change 
due to carbon dioxide (CC/Hr.) 


Reference Times 
level (meters) 0635 0835 1035 1235 1435 1635 


8 -0.10 0.17 0.26 0.49 0.31 0.13 
17 -0.02 0.06 0.15 0.26 0.25 0.08 
35 0 0.03 0.07 0.09 0.10 0.01 
51 0 0.02 0.05 0.06 0.06 0 

100 0 0.02 0.02 0.02 0.02 0 
165 0 0.01 0.01 0.01 0.02 0.01 
240 -0.06 0 0.01 0 -0.01 0 
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4.3. The local rate of change of temperature. 

The local rate of change of temperature was obtained from the MIT 
data at 4.0, 8.0, and 16.0 meters and from the aerograph temperature 
data at 35, 51, 100, 165, and 240 meters. Since the observation times 
were centered bi-hourly at 35 minutes after the even hour, E was taken 
over a two hour period and centered at a time 35 minutes after the odd 
hour. For example, the local change from 0635-0835 CST was centered 
at 0735 CST. The aerograph temperature changes were obtained in a 
like fashion except that some difficulty was encountered in the use of 
this data. A plot of L-20 data with height yields two curves since data 
is taken on both ascent and descent of the aircraft. As an illustration, 
fig. 3 gives the L-20 temperatures for the 1235 sounding. Where data 
for the given reference level were missing smooth-curve interpolation 
was used. This would give an estimated Ve for a two-hour period either 
from the ascent or descent curves or both. Where doubt existed as to the 
representative values of delie. whether to use the ascent or descent 
curves), either an average value from both ascent and descent curves, or 
a subjective choice of the most reasonable sounding was used. Since all 
periods were taken as two hour periods, an error in ۵ ¿may possibly 
exist due to a displacement in time away from the supposed observation 
times. This type of error is estimated to be no more than 5% however, and 
should be 5. in nature. Table 3 presents the observed (or estimated) 


values of It : 
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Table 3. Average local rate of change of temperature C c/nBr) 


Reference Times 
level (meters) 0735 0935 1135 1335 1535 


8 3.24 — 9:27 1.61 1.07 0.12 
17 2055 222 460065 114 Mi 
35 2.08 267 1.41 0.80 0.45 
51 IO ІШ: 22 

100 0.51 2.61 0.87 0.80 0.40 
65 -0.20 2.15 0.93 0.90 0.37 
240 0.15 1.42 1.05 0.97 0.22 





4.4. Advective temperature changes. 

Generally speaking, in computing Ku the assumption will be made 
that advection is zero. However, it may be shown from the quoted geo- 
strophic winds /10/ that this assumption is not quite valid for the period 
in question. Haltiner and Martin Ay have shown that the layer mean 


geostrophic advection is given by: 


с MEME E DS. (18) 
д p iv E» NA 


where Ма is the layer mean geostrophic wind; ШІ, is the layer mean virtual 
temperature; and ДА is the cross-hatched area of the geostrophic wind 
hodograph in fig. 4. 

Now computing the geostrophic mean advection for the layer 950-900 
mb for 0900 CST on 7 September we get 0.24°C/Hr. Since the actual wind 
component perpendicular to the mean isotherms is about 17% of the mean 
geostrophic, the замес й term of 0.249 C/Hr. will be reduced to 
0.0419C/Hr. Since the wind speeds increase in the afternoon it is pro- 
bable that a somewhat larger advection will take place at this time. Un- 
fortunately no other useful geostrophic components are given in the data 
Other than for 0900. This consideration of advection, therefore, requiring 


more extensive investigation will not be pursued further in this paper 


NS for the 0935 sounding where the advection computed for 0900 
will be used in the determination of a К for comparison with a Key 
computed without the advection consideration. 
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4.5. The coefficient of eddy conductivity. 

Having computed (Do, it will be necessary to obtain (ат, Ya) 
Values of (AT 4 Xi) were obtained by adding 0.00989 C/m to aT 
To obtain AT above 17 meters, plots of aerograph temperature with 
height were made on semi-log paper (see fig. 3). 

Similarly for values up to 17 meters, plots of MIT micrometerological 
temperatures with heights were made on linear graph paper and a smooth 
curve drawn through the points. Values of At were then obtained from 
the applicable curve by measuring the slope of the chord centered on the 


l and below, one 


level in question. For vertical gradients at 17 meters 
meter intervals were taken. For the levels 35 meters to 165 meters, 10- 
meter intervals were used, while 20-meter intervals were used above 


165 meters. 


Returning to eq. (6) and omitting advection, one obtains: 


(ar (т ) ] (ar. ) 
К -= [您 )- Be jean | Ni + K< да И ۱ (20) 
_ - —_—_ — A ----- سس‎ 
АТ 4. ) 
2-5). 
Because of the fact that AT was centered after the odd hours and 
Де 


the radiative temperature change was centered after the even hours, it 
was decided to compute a Кн centered after the odd hour by doing the 
following: 

(1) Averaging (т) for two adjacent after "the even hour" 


RAD 


soundings, giving a (£) centered after "the odd hour". 
| | Е / АР 


l.The 16 meter change was assumed to apply at 17 meters. 
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(2) Repeat step (1) again for ۵ ب‎ У, їо obtain (27 um ۳ 


values of which are listed in table 4. 


The previous steps yield the final recursion formula 


Күз L& - iue hte, n (ам) (21) 
un 1 X) 
where Ki j K, are also centered after "the odd hours". 
Table 4. Average lapse rate of potential temperature, (°C/meter) : 
Reference lime 
levels (meters) 0735 0935 1135 1335 1535 
4 0.0088 -0.1011 -0.1301 -0.1292 -0.0992 
8 0.0173 OR 020 7 -0.0402 -0.0482 -0.0287 
17 0.0218 O 20702572 205097892 -0.0082 
35 0.0138 -0.0157 -0.0047 0.0072 -07 0092 
eu 0.0183 -0.0057 -0.0042 -0.0052 -0.0062 
100 0.0178 020023 -0.0042 = ОШО -0.0042 
165 0.0208 0.0073 -0.0027 0.0008 0.0003 
240 0.0270 0.0143 0.0014 0.0003 0.0018 


As an illustration of the form of this equation for obtaining Ko at 17 


meters, Кі; , the equation is: 








AT [NT Es 
= | At Lond E Jennys А? до + Ка „+ За), (21) 
AT Yi) 
Az 17 


with A2=17 meters - 4 meters = 13 meters. 
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Now a starting value of к. at four meters, with temperature changes 


at eight meters, yields a value of КН at 17 meters. Next to obtain the 35- 


meter Кн, а starting value of К а: eight meters with temperature changes 


H 
at 17 meters, is required. After this the computed value of Ko at 17 meters 
is used to get a value at 51 meters, and the scheme is continued in leapfrog 
fashion until all values have been obtained. 

As previously noted in section 1, table 5.a of /10/ gives values of 


2 


min | ). These values were obtained from the correlation of 


Q (mcal cm 
rapid fluctuations of temperature and vertical motion using one-second 
samplings of fast response probes together with the Reynold's formula: 
AAA 

As a general principle, following Priestly ІП the turbulent flux 
of heat may be considered constant with elevation in shallow layers above 
about one meter. The reported values, table 5.a of /10/, seem to sub- 
stantiate this result up to 12 meters . Now values of Q given in table 5.2a 
of A0 / are for 3.0, 6.0 and 12.0 meters. Mum these to get represen- 
tative values of Q in the layer 3 - 12 meters gives starting vales © RH 4 


апа К obtained by using eq. (1). 


H,3 
One further complication exists, however, in the 1235 and 1435 
soundings for which no MIT values of Q are given. Here it was decided 
to estimate values of Q from a comparison of the MIT values of Q and 


theoretical model values of Qg at the earth-air interface. Two such 


theoretical modells: were available with values for the 1235 and 1435 


1 E 
For a discussion of these theoretical models see Chapter 7.3 of /10/ 
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*5 ; one due to Suomi A0/ and one due to Lettau A0 / 

Considering the entire seventh observation period,Lettau's values 
of Q, were found to correlate more highly with the reported values of 
turbulent heat transfer Q than did the Qo values of Suomi. On this basis, 
a ratio of the.average observed Q to the Lettau theoretical value Qg for 
the 1035, 1235, and 1435 observations for all observation periods at 
O'Neill, Nebraska was calculated. 1035 values were included to give 
enough values on which to base this ratio. The ratio was then applied 

to the Q, of the Lettau theory for 1235 and 1435 so that Ол апа Оз for 
these times could be obtained. Computed values of ie are given in 


table 5. 
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S. Results. 
5.1. Ku as a function of height. 


Considering the values of КЕ given in table 5 up to the point where 


there is a sharp drop off in values with respect to height , the computed 


m 


values of Кү show a strong linear trend when plotted on log-log paper. 

To substantiate this linear trend,correlation coefficients of log 
Кн and log Z were computed and regression lines of best fit were obtained 
and plotted for the data for each sounding (see figs. 5-10). 

Now a straight line plot of Ko with 2 on log-log paper yields K 


H 


as a function of exponential height according to the general equation: 
“ут 
Ky = c3 (22) 
Ki “| 
Computed equations for ра as an exponential function of height are 


presented in table 6. 


Table 6. Кн as an exponential function of height. 


т” 
44 


Range of Heights on 


Time (CST) Equation which the equation is based 
(meters) 
0735 Ky = 2602: 4 - 320 
0935 К = a 4 - 0 
09352 К. = 57997 2: 4° | 4 - 320 
1135 Қ, -226e( z < 4 - 165 
1335 Ki, = 22 ии 4 
1535 Ka = roon 4 - 100 


lNegative values, and values above 165 meters for the 1135 sounding, 
were excluded. 


Advection neglected in the computations except for the second result at 
0935, for which the 0900 computed advection was assumed to apply at 0935. 


31 









| 1 ۳ 


iu ро WO р .. a" 


ат (осон хс жі дас сой юз меду T аам во wn Они, 





Series mt кеші ^v Knee ‘`` ۱1 ve yn? Pore Lj] 4 е! 





OER poL-eof re: bee rye Loses ead ees £m „Я $ 
мо аме FI tenons ges. LUV 0095164 [Sag OF 





MEFS a aw ۱ ۱ 00 4 0 Ri) Holes ecu LE DURO e Sal rae I 1 


— 
¡fol E сай SD ahaa a maes mal ek al ات‎ 


SA asqeu mu! тте М ы” E MES mH Mr now 
асите ләҹә эму тв онар Agad сеге і «Шана e 
| | ۲ 业 
ч ЗА x 













DP AA ARA Nr ж. зЗ ۱3۲ 2۰ penam 






n akaw 11) ж 
A ۷ Зі d ^. ۱ VS EPA metes m кə A 47 EU 
m 


ETT ume zu CEREUS ПЭШ жат 










"Pete ip. жа 
kz s. S 
"m кв.- , 3 "ын 
ы 
پت‎ 













"£ w 5 





5 









EE 7 РИС 
жм жы” С 
и. фр») TP А 
em ТҮ ۱ A 


А J = 
к 
i ۱ 


EJ 
5 
B 


а ۱ P | 


Bernier ERES НЫ SEP ASE TPH CHER FE (HEE B > 










































































| ПІШЕН В ии ү 1 “E 
EL | HERE ШШШ a pee 
ol 2 HERE ВЕБ, | ЕЕ 

Dn 1-1 ۱ oo 
il EN: ІШІ ا‎ | 1141 ЕЕ 
ІШ da ү | ОШО ۳ 
ш 
144,2: 
ИНЕШ ШЕН ГА NUES == ES ИШ: - 一 一 
ПІ 1. ІШІ EEE 
ШШ | un ЖЕНЕ | 


ИХ 
El Е qut cR bo c 
ВИЕ EN = 


HERI 
I | 


сесия 
НҮ 
IH: 





БЕ 


Í 
| 


3 


sss,‏ و 














LS. 
za 
up 
DS 
2 








ЗЕ 





ТЕ: 
MEER mm dE 
TT mm КЕНЧ ЕЕ 
108 НІНЕ ها‎ 
E » = 








j ed ЕЕЕЕЕЕЕ 
i But — 
pr 


yla 
ey: 
ін m c 
nd = 
: 
E 
13 т 





Ki as a function of wien i, the ЗЕ S 


— 


ее оь 





ЫЫ = il 



































ii | БЕІН! a ` 

ШЇЇ d un E ae Т 

|| PREECE EEHEG 1. E E Es 
"m | m = 1: 
ы МЕНЕЕ 1 in PE E | ШШШ Г 1 ur р 
IE ایا‎ 
ی‎ И: ШШШ 
Er ml з 
. шш O 
ЕЕ а PA 2 a: З 
р 
| Ve КІ Е ҮШ D EE 
ІШЕ HI HEHEH p MEE Е == | 11 H c ШЕН НЕЕ НЕЕ ы 
ПЕШКЕ НЕ d | НЕГІ | i 
(ШИШ [| ІШЕ tn i "T | 5 
۱ | и | . M | M da 二 Í. 

x 

















1. 


ЧИЕН 


1: x 

IE 

` = I PE 
rH a 


|] "E ү a НЕЕ 


= о a 






pe 

二 全 一 s 

ф ши >< Ies 
2 

E 


| { 
i 1 pt 
۳ IE 
Tatu 
I 
Tib 


КА 
г 
一 
ШЕ 
иы 
ELS О 
ШИ ESEEEES 
В К ERES 
BAL ба 
Т 
Пт — —— 


НЕЕ 
РЕНО 

НЕ - 
E аара аарц 
t ill 6 

55 1 2۵ اش‎ а 





Bas 





- D HH | 
{| | EF |. 2 d oO | 
П Tr uc T. aa ш 
| S EI 51. V T 

ШШ 


|. A n ШІШІ illi pd o) i i 


sl ІШІ 


оо 














HH t Hr и ШЕ ЕЕЕ е ТН ы ed 
| ЕЕ 
ll 


E in => 
Hi Е ll ا‎ 141 





+A 


г 
об 二 -二 一 
TE ЕН 


EE 


ИЕ 










































5 t Е: 
A i | dd Faris ЕЕ EERE Ed үр ШЕ 7 
i s EE sP s НІНЕН В 
N о E 
| Е e 
1 Г "E 12 0: 
Boo o ol 
EN EDA Pe HE tum kR. Е FIS d 
| ү 4. ла E ii il 
| ү he e 
17 ІН | E 0 1 由 ШЕ E E a ۳ z F 
| a | ра | 11. Е: 
D a ов 
2. pr ess ا‎ 
114 12. EI EL us 
| 1 | د‎ 
| || ШІ ШЕШІ i ЧИШТЕ! ІШІ ІШІНІНШІ | a Ще 
D. T DN | DELLI LIE a 
D EER ا‎ Т 0 1 
۳ سس‎ : 
B RU EIC Wee XB Hl ІШІ WE E le 
с: ۳ - ۳ 
B | 38 01 

BUE ү ae EN ШЕЕ: de 


AE see - | 1. 
14-11 | 
ІІ RE E 














SITE ҮШҮҮ 
ШП Т s Е rud 


Enn Т | 1225 ADD HIHI 


1 


t 


— =e 
— ~ 十 -一 
To > 


4 
HII 


| 


= 







































i | ji OT | ii 
Ee | | cu E 1. | | || 
BERTA INS EE | ШШ я. ШЕШ d id 
N 
x i | | Ош | ВЕ. 
o ll 
|| || | | Ш ІШІ 11710 || = | | 
Е | Е 1 IE min 
E ТОП HERE R TE àv ТЕ EE 
ІШЕ 111. d ПЕЕ | | a | 
141071 1524 
ІШЕ 53 ШШШ ШО 11 i ШЕПТЕ 1 | 1. ІШ (ТІГЕ He 
271. E 
171... 
ll o 
1 | x | | ا‎ 
ҮШҮ НЫНЕ 
11111 ۱ | | | || | ИН a 
ІШІН ШТ ТЕ IE 
EHE р I 1 LEIDEN | | 
luas Il l) 2. Ша. 
HA 
۱ ЧЕНЕСЕ ШШЕН. ud E ЕЕ. ШІ uc ШЕ E ЕЕ "E 
. p 111. о ان‎ ге. la ju 
Dedi -8 72 ШІ Н Е ШЕШ ІШІ Bis 
ШЕ с | = B m ni ү im "LS |. ШШШ HE p EUH 2H 
Е (n. a m „+ Е ШІ СИ d 
НЕЕ | الل ا ا‎ 





r 


Kg as a function of height for the 1135 sounding. 


B WE 


Fig. 











š TN ТП ЕЕ: E‏ و 
cs pug PN pen eres a Fer Eon ШЕ qe‏ - 
Wu. ы p. 7 ма Eu‏ 
PE 1‏ 






















o 


а = E Е ўа = 0 === = = ERE 
= بت‎ 


iE Давай нат ва 








- 1-2 > 
= А zr : 
F 2 e == 


тер? Е 
= Е. ~ -~ 一 


A A 一 























EL 
ap 
= 
Lala] 
== 


| ТІ E 
1. 
»" ; | | ۹ 


Е 


de ШІН DEN 
5 | ШШШ i iB pue = 


il | 4. ШІ i || 
17 ub 11 

| Ht ШІ | НІНІ m 
| | i ИШИН | 


| | 
EORR IRI RES 

ДЕ ОАЕ ВЕНЕРА ЕНЕ E 
ep 1 ШЕННЕН 
E H Ци ІШ ІШІНЕ “| 
E ЕНТ! ШІ ЕНЕГЕ. 
d NI || EE HII EH TS = ШЕ 


|. | p E 
17 2| 


1 1 El: 
a ү || ПТ | ШІН ІШІН: | 
Hi 
Mix: 
[IBS 
3 E 
pr E i Uu H = 
| ND | — 4. 
р | ЕНЕВ نی‎ 


| г Т! a 
qx m Ug 


HI ao el; 
Er lie dd Rr 111. QUNM e 
Н. Í DN Там i i ТЕ ad III | 

M trin 


ІП ЕТТІ 


1 
ii | ІШІН ІШІ НЫ ч 1 
Paese ЕЕ 


CE 





=== 4 
ЕЕ 
BESE PCIE 
o 
i. 
Е ЕЕ 
a EEE КЧ 


<< 228 


E 


SSS gE 
ЕЕ Е er 


[pd on pa دس‎ 








=— ey 
=o 
eer 





| | Ін 
HH 
|! | na 
ШЕ | 
= ПШ 
E 4 ۳ 
ux i isa 








ЕЕ 


- sounding 


4 


tii 
H: 


E 


us 


4 


ج 
к» |‏ 


ШТ 
Ha 


ale 
i 


24—09 


| 









ES 
A 
жол: 


pred 


DI DTI — 
| 合计 天 只 站 加 الب لس‎ 
GIZ Em = SSES eL Ep 
ا که‎ 
Кыскы е аа 


== 






تسس | 


НЕТ 


реда 


ЕЕ 
РЕ == == ЕЕЕ 


ЕЕЕ ЕЕЕ НЕТ 
Е Ц 
[ТЇ БЕШ ШЕ ШИ БИ ШИ ШШ ЕШ о man 


ins 
Mu 
ІНІ 


= 
c ben: 
ШЕЕ Г 


























аве 

e 1 

Soe BESE 2222 
== F 

EE rri ПІН 
Sssass 8 Г: 

FFF С 

КЕ е а 1 pu 

peu] 

1 5 8 8 7 11 





























ШШ ШІ ІШІ 
ШІ ІІ 11: 
i ІШІ ІІ | QE 
ЕЕ ШЕШШ: | ШЕ 
"T mI | 


LLL 


11! 


ONU CEDE 
ш к= = 


> as a function of height for the 1 





F 
IEA 
HHH 
E 
HH БЕ 


št 

dH 

н о В 
H ү HE 
uu us i. 


т 
ЕҢ 
ИТ 

< s 
E rtr 
------ 
BE E Û RY. 

ЕЕ‏ ز225 


FIS 10. 


a || E 3H 


ЕЕЕ Е 


ys 
ш 





- It ıs to be noted that where advection was considered (0935 case) 
the result was not greatly different from that in which this parameter was 
neglected. The dıfference, however, indicates that the ре values, 
with advection included, were slightly greater in the lower portion of the 
layer, and smaller in the upper portion. 

It should also be pointed out that in most cases, the values of En 
computed for elevations above 165 meters showed a decrease below the 
value at the previous level. In several cases, the values of Ке уеге 
negative. This was most likely due їо observational errors in (TA. +54 ) j 
when 01/۶ was very close to the adiabatic lapse, a condition which 


occurred in the neighborhood of 165 meters for the midday and afternoon 


ЕМТІЗ , 
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ша. Кн as a function of stability. 

It would have been appropriate to compare values of Kas and m 
to a stability parameter such as the Richardson number. Since appropriate 
wind data for such a comparison were not available , it was more feasible 
to correlate ка and m with a significant part of the 222 of the in- 
tegrated Richardson тин (NC . This was done and the results are tabu- 


lated below. 


Table 7. Linear correlation coefficients of AO with K and AO with m 
at indicated heights. H 


E (meters) AG computed Correlation Correlation 
over layers(meters) coefficients of coefficients of 
Ku with AO m with AC 
4 2-4 0.26 -0.73 
Ө. 2-8 -0.77 -0.71 
17 2-17 -0.96 -0.68 
39 2-35 -0.88 -0.84 
51 2-51 O DS 
100 2-100 -ü on OPE TES 


Considering values of AQ versus Ky at heights from four to 100 
meters for all soundings (sample size 30) a correlation coefficient of 


-0.43 is obtained. í 


Values obtained from equations given in table 6. 
Неге ЛӨ = Ө(2)- Ө(2) 


where Z is variable and ©(2) is the value of G at 2 = 2 meters. 
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Ga, Conclusions. 


Although negative values of K.. are somewhat disconcerting, other 


H 
investigators have been confronted with this anomaly at much lower levels 
/ 6 / As mentioned earlier, factors possibly contributing to these negative 
values are chiefly non-representative lapse rates due to errors in the 
measurement of temperature, especially as the lapse rate of temperature 
approaches the dry adiabatic. Still another possibility arises when one 
considers the theory of Priestly and Swinbank concerning forced versus 
free convection (see for example p. 261 of Р à 

From an inspection of the computed power-law profiles for кё listeq 
in table 6, it appears that m increases inversely with stability, the maxi- 
mum value of m occuring in the early afternoon. Although it is doubtful 
whether the correlations listed in table 6 can be proved statistically sig- 
nificant on the basis of such a small sample size (five), nevertheless it 
is felt that a strong physical significance is inherent in the repetative 


nature of the results, and that K.. has a strong inverse dependence on 


Н 
stability. This statement is given support by the results of a correlation 
of the AS parameter with all values of Кн represented in table 7 (sample 
size 30) where a statistically significant! correlation coefficient of -0.43 
is obtained. Further it is apparent, from an inspection of table 7 that this 
correlation (-0.43) could be further improved by the elimination from con- 


sideration of values of КЫ at four meters where a positive correlation of 


K., . with AO is abtained. Considering the positive correlation of К 


Hw H3 


leat the five percent level of belief. 
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with ^.Q and the negative correlation of ки 0 with AQ , it seems likely 
that the method of computing КН at four and eight meters may have affected 
the overall correlation. It must be remembered that the assumption made 

in this connection was that Отокв is independent of elevation from one 

to 12 meters. At higher levels, the computational procedure gave con- 
siderably larger Кн-уашез which then varied consistently (if not "sighi- 
ficantly") with stability. 

In summary , daytime-values of Kr; obtained under the conditions of 
the Great Plains Turbulence Field Program ANA tend to increase with іп- 
creasing lapse rate, and with increasing height up to at least 100 meters. 
Above some such level КН seems to decrease. It is possible that a 
specific function Ка = £ (Ri) may exist, but the data was not sufficiently 


detailed to establish such a result at this time. 
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